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Abstract -A measure of quantum correlation defined from an information-theoretic perspective, 
namely, quantum discord, is applied to study the time-evolved nonequilibrium state of the infinite 
anisotropic quantum XY spin chain in a transverse time-dependent field. In particular, we probe 
whether the collapse and revival of nearest-neighbor entanglement of the state seen with a varying 
initial applied field strength, at a fixed evolution time, may be predicted from the behavior of 
the quantum correlation measure. For this quantum many-body system, realizable with currently 
available technology, we find that the revival of entanglement of the evolved state happens if there 
is an increase in quantum discord in the vicinity of entanglement collapse. 
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Introduction. — There have been extensive studies 
in many-body and quantum information theory (QIT), at- 
tempting to quantify the key aspects of quantum correla- 
tions present between the parts of a system. The foremost 
among them is on the entanglement-separability distinc- 
tion and its subsequent quantification by using measures 
such as entanglement of formation, distillable entangle- 
ment, and relative entropy of entanglement (see [1] for a 
review). Entanglement has been extensively used for in- 
dicating quantum criticality in phase transitions [2] and 
in numerous QIT applications such as quantum telepor- 
tation [3], quantum dense coding [4], and quantum key 
distribution [•">]. However, there is no quantitative agree- 
ment between the various measures of entanglement [1]. 
Even qualitative differences appear: Distillable entangle- 
ment vanishes for a certain class of entangled states, called 
bound entangled states [6], although other entanglement 
measures produce nonzero values. 

Moreover, there exist quantum correlations which ap- 
pear even when entanglement is absent. Several phenom- 
ena have been discovered which produce nonclassical re- 
sults with no shared entanglement identified in the system. 
These correlations feature in important aspects of QIT, 
such as "quantum non-locality without entanglement" [7] 
(cf. [8]) and "quantum data hiding in separable states" [9]. 
A natural question then arises: What form of quantum 
correlations is responsible for such nonclassical (nonlocal) 



behavior even in the absence of entanglement? 

Early attempts to explore the concept of nonclassi- 
cality and correlations from a perspective that is differ- 
ent from the entanglement-separability paradigm include 
those defined in the language of quantum optics [10], 
and the literature on Bell inequalities [11]. Recently, 
information-theoretic and thermodynamic concepts have 
been used to define quantum correlations independent of 
the entanglement-separability criterion - the measures of 
quantum discord (QD) [12, 13] and quantum work-deficit 
(QWD) [14] quantify quantum correlations by attempting 
to quantize expressions for correlations existing in classi- 
cal information theory. The mapping from a classical to a 
quantum system often introduces a lack of concord, mainly 
due to the non-commutativity of quantum operators. QD 
quantifies quantum correlations by using the difference in 
the quantum expressions corresponding to two equivalent 
definitions of classical mutual information [15]. QWD, 
on the other hand, does it by using the difference in the 
amount of negentropy ( "work" ) extractable by global and 
local heat engines [16, 17]. The application of such mea- 
sures of quantum correlations in many-body systems may 
reveal new phenomena which cannot be detected by entan- 
glement. QD has already been applied for studying static 
properties, such as properties of the ground state and the 
equilibrium states, of quantum spin systems [18]. There is 
currently no concrete operational relation between entan- 
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glement and other measures of quantum correlations such 
as QD or QWD. 

In order to probe this interrelation, we focus here 
on a known, exactly solvable model, that of an infinite 
anisotropic XY (quantum) spin chain in a transverse field. 
The nearest-neighbor entanglement in the time-evolved 
state of the chain, at a given fixed time, exhibits criti- 
cal behavior - a dynamical phase transition (DPT), con- 
trolled by the initial value a of the transverse field [19] . We 
study the dynamics of quantum correlations, as quantified 
by QD between two neighboring spins, and show that, for 
this particular system, the observed collapse of entangle- 
ment and higher-field revival (with changing initial value a 
of the transverse field) can be predicted by the behavior of 
QD - entanglement revives after collapse, if QD increases 
in the vicinity of the region where entanglement collapse 
occurs . Hence, an increasing QD at the point of entan- 
glement collapse (i.e., at the corresponding time and field 
strength) can be used as an indicator for the entanglement 
revival. We further show that this behavior is potentially 
generic in that it does not depend on the specific measure 
of quantum correlation used - QWD also shows a simi- 
lar predictive capacity, although further work is needed to 
extract a straightforward criterion [20]. 

Measures of quantum correlation. For com- 
pleteness, we first define the correlation measures. 

Quantum discord. The total correlation in any bipar- 
tite quantum system, pab shared between two parties A 
and B, can be measured by using quantum mutual infor- 
mation [21] (see also [17,22]): 

I {pab) = S{p A ) + S(p B ) ~ S{pab), (1) 

where S{p) = — tr(p log 2 p) is the von Neumann entropy of 
the quantum state p, and pa and ps are the local density 
matrices of pab ■ 

Classical mutual information can be similarly defined 
for a joint probability distribution {pij} as I{{pij}) = 
H{{ Pi .}) + H{{p. 3 }) - H{{ Pl] }), where H{{q 3 }) = 
~ Ylj 1j 1°S2 1j i s the Shannon entropy of the probabil- 
ity distribution {qj}, and {pi.}, {p.j} are the marginals of 
{pij}- There is an equivalent classical expression for mu- 
tual information using the concept of conditional entropy: 

H{{ Pij }) = ff({pj}) + fr({ P4 | J }) = m Pi .}) + H{{ Pjli }), 

where {pi\j} and {Pj\i} are the conditional probability dis- 
tributions. 

For our bipartite quantum system, measuring on the 
subsystem B using the set of projectors {B{\ [where 
B[Bj = SijBi, J^i^i = ^B, with Ib being the iden- 
tity operator on the Hilbert space on which ps is de- 
fined], when the two-particle system is in the quantum 
state pab, produces the post-measurement states p\ B = 
j-Ia ® BiplA ® Bi, where 1a is the identity operator on 
subsystem A, andp^ = tvAB(lA®BiplA®Bi). The condi- 
tional quantum states that are produced at A, due to the 

1 The collapse and revival of entanglement occurs with the chang- 
ing initial value a of the transverse field for a fixed evolution time. 



measurement at B, are p A \% = ^rtr B(J-A®Bip\A®Bi) , with 
probability pi . The quantum conditional entropy can then 
be defined as S{pa\b) = mm {Si} J2iPi^{PA\i), and simi- 
larly, S{pb\a)- One can quantize the classical expression 

I{{ Pij }) = H({ Pi .}) - H{{ Pi]j }) (2) 
for mutual information to obtain the following: 

J{pab) = S{p A ) - S{pa\b)- (3) 
QD is then defined as 

Q{pab) = I {pab) - J{pab)- (4) 

QD is positive for all quantum states. 

Quantum work- deficit. The concept of QWD is based 
on the fact that information can be treated as a ther- 
modynamic resource [23]. Given a quantum state p, one 
defines the allowed class of global quantum operations, 
called "closed operations" (CO), as arbitrary sequences 
of the following operations: (Gl) unitary operations and 
(G2) dephasing p by using a set of projectors {Pi}, i.e., 
P -> EiPipPh where P.Pj = %P;, Y,i P i = I with I 
being the identity operator on the Hilbert space 7-L on 
which p is defined. Under this class of operations, it can 
be shown (see [14, 24]) that the number of pure qubits 
that can be extracted from p is Ig{p) — N — S{p), where 
N = log 2 dim%. 

Correspondingly, the allowed class of local operations is 
called "closed local operations and classical communica- 
tion" (CLOCC), and is defined as arbitrary compositions 
of the following operations: (LI) local unitary operations, 
(L2a) local dephasing and (L2b) sending a completely de- 
phased subsystem from one party to another over a noise- 
less quantum channel. Let us consider a bipartite quan- 
tum state p = pab ■ The number of qubits that can be ex- 
tracted from a bipartite quantum state pab under CLOCC 
is 

Il{pab) = N- inf [S{p' A ) + S{p' B )}, (5) 

AeCLOCC 

where p' A = tr B (A(p^ B )), p'b = tr A {A{p A B), and now 
N = log 2 dim Hab , with Hab being the Hilbert space on 
which pab is defined. QWD is then defined as 

A (pab) = Ig{pab) - Il{pab)- (6) 

Entanglement: Logarithmic negativity. There is a 
plethora of entanglement measures that is known in the 
community, and each has its own operational relevance 
and degree of mathematical tractability. We will use the 
logarithmic negativity (LN), which is a very useful com- 
putable measure of entanglement [25]. 

The definition of LN is based on the fact that the nega- 
tivity of the partial transpose of a bipartite quantum state 
is a sufficient condition for the state to be entangled - the 
Peres-Horodecki separability criterion [26]. Moreover, for 
two-qubit systems, which will be our domain of study, the 
condition is necessary and sufficient [27]. 



p-2 



Quantum discord surge heralds entanglement revival 



LN is evaluated by using a quantity called "negativity" , 
defined for the state pab as 



AT(pab) 



Pab 



1 



(7) 



where 

Pab of Pab 



Pab 



is the trace norm of the partial transpose 



From the Pcrcs-Horodccki separability cri- 
terion, the partial transpose p T A A B should be positive for 
all separable states. Hence N(pab) is zero for separable 
states. LN of pab is defined as 



Eat (pab) = log 2 

= log 2 [2A/"(p AB ) + l]. 



Pab 



(8) 



The system: infinite quantum XY spin chain in 
a transverse field. — The infinite anisotropic XY spin 
chain in a transverse field is governed by the Hamiltonian, 



H= JY,[(1+7)S*S? +1 + (1 - 7 )SfS y 



i+li 



(9) 

where the anisotropy 7 is nonzero, and J measures the in- 
teraction strength. Si = (J=x,y,z) are one- half of the 
Pauli spin matrices at the corresponding site. Note that 
7 = corresponds to the XX model while 7 = 1 is for 
the Ising model. Here we consider the models for 7 > 0, 
so that the interaction and the field parts of the Hamilto- 
nian do not commute, whereby the external field can have 
nontrivial effects on the evolution. The transverse field is 
applied in the form of an initial disturbance: 



h(t) 




(10) 



where a^O. 

The above Hamiltonian can be realized in a system 
of cold atoms confined in an optical lattice. The two- 
component Bose-Bose and Fermi-Fermi mixtures, in the 
strong coupling limit with suitable tuning of scattering 
length and additional tunneling in the system can be de- 
scribed by the above Hamiltonian [2,28]. The dynamics 
of the system can be simulated by controlling the system 
parameters and the applied transverse field [29]. 

Suppose that the system starts off from the initial state 
which is a (canonical) equilibrium state at temperature 
T. We are interested in the nearest-neighbor (two-site) 
density matrix of the evolved state at time t, that started 
off from the equilibrium state. In general, a two-qubit 
density matrix is of the form 

E ]=x ,y, z M*(t)(fX* ® I + I ® (j) 

where T^ k (t) are the two-site correlation functions, and 
M 1 * {t) are the magnetizations. Using properties of the 



XY Hamiltonian, some simplifications can be made, and 
the final form of the two-site density matrix is [30-32] 



Pf{t) 



I + M z {t){a z ®I + I®a z ) 
+ T xv (t)(a x ®a y + a y ®a x ) 
+ T*>(t)o*®ot). (11) 



Diagonalizing the Hamiltonian via Jordan- Wigncr and 
Fourier transformations, the correlations and the trans- 
verse magnetization in (11), for an initial temperature T 
= 0, arc found to be [30-32]: T x v{t) = Ty x (t) = S(t), 
T xx {t) = G(-l,t), T yy (t) = G(l,f), and T zz (t) = 
[M z (t)] 2 - G{\,t)G{-l,t) + [S{t)f, where G(R,t) (for 
R = ± 1), S(t) are given by 



G(R,t) = 1 [ d<t> BinUR) sin 0— -j- 
7T J A(a)A 



)A 2 (0) 



x \j 2 sin 2 <j> + (cos 4> — a) 
+ acos0cos[2A(O)£]} 



COS ( 



1 r 1 

/ d<t> COS 6 — r—TT, x 

tt Jo V V A(a)A2(0) 
x ({7 2 sin 2 (f> + (cos <f> — a) cos < 



COS( 



a7 2 sin 2 <^cos[2A(0)t] 



S(t) = 



1 a 



and 



M z (t) = - 



d(j> sin 



1 



, sin[2tA(0)] 
' A(5)A(0) ' 



o "A(5)A 2 (0) 
{cos[2A(0)f]7 2 asin 2 <p] 
cos 4> [(cos (j> — a) cos (j) + r y 2 sin 2 <j>\ } . 



Here A(a;) 



b 2 



+ [x — cos (j)] 2 } 2 , and a 



a/J, t = Jt/h. We will use d and t as the (dimension- 
less) initial field and time parameters, respectively 

Methodology. — We observe the behavior of the sys- 
tem at a time t after it starts from an initial canonical 
equilibrium state at zero temperature. In particular, we 
wish to find the measurement strategy for obtaining the 
optimal QD for the nearest-neighbor (two-qubit) density 
matrix of the evolved state at time t. Since we focus on 
projection- valued measurements, and since the local sys- 
tems are qubits, the measurement will necessarily involve 
projecting onto an orthonormal (two-element) basis of a 
two-dimensional complex Hilbert space. Let that basis be 
given by 

= cos||0) + e^8in||l>, 



*2 



e-^sin 



-|0) + cos-|l), (12) 
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where {|0), |1)} form the computational qubit basis. 

In the case of QWD, the general definition involves 
two-way communication of dephased states, and is as yet 
not computable for arbitrary states. We consider the re- 
stricted case where only one-way communication is al- 
lowed. Again, as the local subsystems are qubits, the most 
general measurement basis for the dephasing will be of the 
form in Eq. (12). 

In our analysis, we consider the state of the system as we 
sweep over the applied initial field strength a, at a fixed 
value of the anisotropy parameter 7 for a fixed time of 
evolution of the zero-temperature equilibrium state. The 
motivation for considering such a state is the fact that 
the anisotropic transverse XY model, at zero temperature, 
undergoes a quantum phase transition (QPT) at a = 1, as 
we sweep over the magnetic field at a fixed 7. An initial 
zero-temperature equilibrium state ensures that thermal 
fluctuations are absent. 

Entanglement versus Quantum Discord. We 

consider the initial states to be equilibruim states at zero 
temperature. The observed DPTs occur for a wide range 
of the asymmetry parameter 7. In Fig. 1, we plot LN of the 
nearest-neighbor state of the evolved state, as a function 
of time t and the initial field strength a. An important 



Entanglement 




Fig. 1: (Color online) Behavior of entanglement (as quanti- 
fied by logarithmic negativity, measured in ebits) with respect 
to the time of evolution t (dimensionless) and the initial field 
strength a (dimensionless). The collapse of entanglement to 
zero are distinctly shown in the regions around t = 1 and t = 4. 
We have taken 7= |. 

observation from the figure is that for times until a little 
before t = 2, and again at around t = 4, the entanglement 
collapses to zero at a certain field value, but revives to 
give nonzero values at higher fields in the dynamic evo- 
lution. At other times, the entanglement recedes to zero, 
and does not become nonzero at higher field values. This 
remarkable behavior of collapse and higher-field revival 
of entanglement in this case takes place with respect to 
the varying initial field applied to the system 2 . For defi- 
niteness, the figure is plotted with the anisotropy 7=5- 

2 A temporal collapse and revival of entanglement can also be ob- 
served in Fig. 1. The entanglement collapses and revives with varying 
time for certain values of the initial field strength (a > 1.2). Our pri- 
mary interests, however, lie at understanding the correlations that 



However, all the results hold irrespective of the value of 7 
chosen in the range (0, 1]. 

How does QD behave in this interesting range of initial 
field strength a and time t? The first observation in our 
computed behavior of QD shown in Fig. 2 is that QD is 
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Fig. 2: (Color online) Behavior of quantum discord (measured 
in bits) as a function of time t (dimensionless) and initial field 
strength a (dimensionless). 7 = |, as before. QD shows a 
different behavior around times t = 1 and t = 4 for which 
entanglement reappearance takes place. This is made clearer 
in Fig. 3. 

nonzero at points where entanglement is zero. We further 
observe that the behavior of QD at points where entan- 
glement revival occurs is markedly different. 

The system thus supports another form of quantum cor- 
relation, even when entanglement goes to zero. This leads 
us to pose the query: Are there more general forms of 
quantum correlations in the system, whose presence at a 
point of entanglement collapse, can predict the revival of 
entanglement for certain times of evolution, and whose 
absence anticipates a non-revival? We give an affirmative 
answer to this question and formulate a relation that as- 
certains the observed behavior for the considered system. 
Consider a time-evolved nonequilibrium bipartite quantum 
state PabO^)' f or an anisotropic XY spin chain in a trans- 
verse field, obtained by time-evolution for a duration t and 
varying with a system parameter a. If for a fixed time t, 
the entanglement E vanishes at a = a c , then 

J)Q(pab) > Q E(pab) > for some \a\ > \a c \. 
oa ~ a c 

(13) 

Here Q stands for quantum discord 3 and the parameter a 
is the initial transverse field. From (13) it is evident that 

arise close to the zero-temperature QPT at a fixed time of evolution 
of the quantum state. Note that the QPT at zero temperature also 
happens at a fixed time, which is t = 0, as we sweep over the a axis. 
The DPT, observed in Ref. [19] , considers the status of the zero-time 
transition at nonzero times, if we still sweep over the a axis. 

3 For the specific model, the violation of the inequality in (13), 
usually results in no revival of entanglement. There may, however, 
exist few intermediary states where entanglement revival occurs in 
absence of distinct positivity and the collapse-revival behavior of 
entanglement remains inconclusive. For positive a such resurgence 
of entanglement without positivity of (13) around a = a c , happens 
in the vicinity of a = 1, and so it is plausible that such exceptional 
cases are related to the zero-temperature QPT in this model [28]. 
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the QD decreases if one considers DPT along the negative 
a axis. This is due to the reflection symmetry and surge 
(dip) in QD is along the positive (negative) direction of the 
considered physical parameter. Hence we consider just the 
positive parameter axis and the surge of QD. 

Already in Fig. 2, we see that QD behaves differently for 
times for which entanglement revival takes place. To get a 
clearer picture of the situation, in Fig. 3, we plot entangle- 
ment and QD for different fixed times t, as functions of the 
initial field strength a. Only for those times t for which 




Fig. 3: (Color online) Quantum discord (continuous, black) and 
entanglement (logarithmic negativity) (beaded, red) for fixed 
times t as functions of the initial field strength 2. Entanglement 
collapses to zero and revives at higher-fields for times before 
t as 2 and at around t = 4. For these and only these times, QD 
increases around the collapse of LN. LN is measured in ebits, 
QD in bits, while the horizontal axes denoting the initial field 
strength a are dimensionless. 7 = ~, as before. 

QD is an increasing function of the field strength a at 
a = a c , a c being the field strength at which entanglement 
vanishes for a given t, there is a revival of entanglement at 
a higher value of a, as stated in (13). In Fig. 4, we plot the 
partial derivative of QD with respect to a, at a = a c . It 
is a pictorial representation of the relation given in (13). 
A similar conjecture can also be obtained by studying the 
correlation dynamics by varying the anisotropy parame- 
ter 7 even though no consistent DPT is observed for such 
evolutions. 

Behavior of Quantum Work-Deficit. We consider a 
second information-theoretic measure of quantum correla- 
tions, QWD. Sections of its plot at different times is given 
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Fig. 4: (Color online) Increase of quantum discord at entangle- 
ment collapse indicates entanglement revival. The continuous 
red curve is of the (first) partial derivative of QD (measured 
in bits) with respect to the initial field strength at a = a c ver- 
sus time t (dimensionless). The black diamonds denote the 
maximum values attained by LN after its collapse (measured 
in ebits) versus t. The continuous curve crosses over zero for 
t < 2, and again at t » 4, which are exactly the times for 
which entanglement revival happens, i.e., where the curve of 
diamonds is nonzero. The inset magnifies the crossing over zero 
around t ~ 4. Thus the maximum entanglement after collapse 
is nonzero only at times where the partial derivative curve is 
positive (since a is positive, the relation in (13) holds). 7 = |, 
as before. 



in Fig. 5, and it is seen that around times t = 1 and t = 4, 
QWD reaches high values for moderately large a. Hence 
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Fig. 5: (Color online) Quantum work-deficit (measured in 
qubits) is plotted against the initial field strength a (dimension- 
less) for different times t (dimensionless). Around times t = 1 
and t = 4, for which entanglement reappearance is observed, 
QWD rises to high values for moderately large a. Again, 7 = 5. 



even though no working relation, unlike for QD, exists for 
QWD, it carries a predictive capacity, which will be ex- 
plored in a future work [20]. 

Discussions. — We have shown that a measure of 
quantum correlations known as QD, that goes beyond the 
standard entanglement-separability paradigm, can be used 
to encapsulate a criterion that tells us when entanglement 
will revive after collapse in an infinite anisotropic quantum 
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XY spin chain. Indeed, we have shown that an increase of 
QD around entanglement collapse indicates a revival of en- 
tanglement at higher fields. This result could potentially 
be useful in applications of entanglement as a resource in 
quantum communication and other quantum information 
tasks to indicate which parameter zones of a given physical 
system contain entangled states, and which do not. 

We would like to test whether the findings are generic, 
in the sense that other measures of quantum correlations 
could show the same behavior, by considering a second 
information-theoretic measure of quantum correlations, 
QWD. We have seen that even though no working relation, 
unlike for QD, exists for QWD, it carries a predictive ca- 
pacity which could be a generic feature of other measures 
of quantum correlations defined beyond the entanglement- 
separability criteria. 

Studies in information-theoretic measures of quantum 
correlations, such as QD and QWD, have revealed that 
such measures give a fine-grained picture of quantum 
states of distributed systems in comparison to that pro- 
vided by entanglement. Our study indicates that such a 
fine-grained picture can show the underlying reason for 
the dynamics seen for entanglement in quantum many- 
body systems. 
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